arXiv:1505.03400v2 [quant-ph] 24Jul2015 


Tunneling time in attosecond experiments and time-energy 

uncertainty relation 

Ossama Kullie 

Theoretical Physics, Institute for Physics, 

Department of Mathematics and Natural Science, University of Kassel, Germany* 

In this work we present a theoretical model supported with a physical reasoning leading to a relation which 
performs an excellent estimation for the tunneling time in attosecond and strong field experiments, where 
we address the important case of the He-atom [1, 2]. Our tunneling time estimation is found by utilizing the 
time-energy uncertainty relation and represents a quantum clock. The tunneling time is also featured as the 
time of passage through the barrier similarly to the Einstein’s photon box Gedanken experiment. Our work 
tackles an important study case for the theory of time in quantum mechanics, and is very promising for the 
search for a (general) time operator in quantum mechanics. The work can be seen as a new fundamental step 
in dealing with the tunneling time in strong field and ultra-fast science, and is appealing for more elaborate 
treatments using quantum wave packet dynamics and especially for complex atoms and molecules. 

Keywords: Tunneling time in strong field and ultra-fast science, time measurement in attosecond exper¬ 
iments, quantum clock, time-energy uncertainty relation, time and time-operator in quantum mechanics, 
photon box Gedanken experiment. 


Electronic mail: kullie@uni-kassel.de 



2 


I. INTRODUCTION 

A comprehensive theory of time measurement in quantum mechanics is missing to date [3] (chap. 3). Often is said that 
time plays a role essentially different from the role of the position in quantum mechanics. In contrast Hilgevoord [4] argued that 
there is nothing in the formalism of the quantum mechanics that forces us to treat time and position differently. Observables 
as position, velocity, etc. both in classical mechanics as well as in quantum mechanics, are relative observables, and one never 
measures the absolute position of a particle, but the distance in between the particle and some other object [4, 5]. Indeed there 
is many attempts to consider time as a dynamical intrinsic, or an observable time called event time. Hilgevoord concluded [4] 
that when looking to a time operator a distinction must be made between universal time coordinate t, a c-number like a space 
coordinate, and the dynamical time variable of a physical system suited in space-time, i.e. clocks. Busch [6, 7] argued that 
the conundrum of the time-energy uncertainty relation (TEUR) in quantum mechanics is related in first place to the fact that 
the time is identified as a parameter in Schrodinger equation (SEQ). He classified three types of time in quantum mechanics: 
external time (parametric or laboratory time), intrinsic or dynamical time and observable time. External time are carried out 
with clocks that are not dynamically connected with the object studied in the experiment, and usually called parametric time. 
The intrinsic or dynamical time is measured in term of the physical system undergoing dynamically a change, where every 
dynamical variable marks the passage of time, we will see this is important for our time invention where the energy serves as 
the dynamical variable in question, which enables a quantitative measure for the length of the time interval of the tunneling 
or the tunneling time (T-time) in strong field experiments. The third type according to Busch is the observable time or event 
time, for example the time of arrival of the decay products at a detector. 

In the history of the quantum mechanics, the earliest attempt, which causes one of the most impressive debates, is the 
Einstein’s photon box Gedanken experiment [8] or the Bohr-Einstein weighing photon box Gedanken experiment [5] (and the 
references therein). A photon is allowed to escape from a box through a hole, which is closed and opened temporarily by a 
shutter, the period of time is determined by a clock, which is part of the box system, that means the time is intrinsic and 
dynamically connected with the system under consideration. The total mass of the box before and after a photon passes is 
measured. Bohr showed that the process of weighting introduces a quantum uncertainty (in the gravitational field) leading 
to an uncertainty in time r, which is the time needed to pass out of the box that usually called the time of passage [6, 7], in 
accordance with the TEUR, eq (1) below. Aharonov and Rezinik [5] offer a similar interpretation, that the weighing leads, 
due to the back reaction of the system underlying a perturbation (energy measurement), to an uncertainty in the time of the 
internal clock relative to the external time [5]. Hence for quantum systems it is important to observe the time from within the 
system or using an internal clock. Busch [6] presented an argument which makes no assumptions concerning the method of 
measurement, and simply based only on a version of quantum clock uncertainty relation as follows, if the energy of the escaping 
photons is determined with an accuracy 5E from the difference of energy before and after the opening period of the shutter, 
then these energies must be defined within an uncertainty SE, i.e the box energy uncertainty AE must satisfy AE < SE. Then 
the clock uncertainty allows to conclude that the box system needs at least a time to = in order to evolve from the initial 
state “shutter closed” to the orthogonal state “shutter open”. Accordingly, the time interval within which a photon can pass 
the shutter is indeterminate by an amount AT = to- This leads to Bohr’s TEUR AT5E ^ h [3] (chap. 3). 

In this work we use similar ideas, where the T-times under the barrier (denoted Tj, is suggested to be similar to the time 
of passage through the barrier (and escaping at the exit of the barrier), and the (quantum) particle (an electron) undergoes 
this process spends a time that is the time needed from the moment of entering the barrier to the moment of escaping from 
under the barrier in the tunneling direction. In addition we suggest a time interval needed to reach the entrance of the barrier 
(denoted Trj, ■), after it is shaken off by the laser field at its initial position Xi. ^ is similar to the traversal time used in context 
of the tunneling approaches [9, 10] or the Feynman path integral (FPI) approach [11-13] (and [3] chap. 7). But in contrast 
we do not make any assumption about paths inside the barrier, where as well-known, the FPI approach is based on all paths 
starting at the entrance of the barrier at t = 0 and end at the exit of the barrier at time t = r, which defines a time duration r. 
A second type of T-time that we invent is what we call the symmetrical T-time or the total T-time (denoted sym)- 
see that can be easily calculated from the symmetry property of the T-time but then later we found ~ i that 

is the time is accounted from the moment of starting the interaction process, where the electron getting a shake-off, responses 
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and jumping up to the tunneling “entrance“ point taking the (opposite) orientation of the field, moving then under the barrier 
(the tunneling) to the ’’exit” point and escapes the barrier to the continuum. The key issue of the present work is (in the words 
of Busch [3]) a case study, the T-time in attosecond experiment and ultra-fast science derived by utilizing the time-energy 
uncertainty relation (TEUR): 

AT- AE> - (1) 

- 2 ^ ^ 

In sec II we present our theoretical model, in sec III we offer a convincing physical reasoning for our theoretical model, in IV 
we discuss our result with a comparison to the experiment and finally we give a conclusion to our work. 


II. THEORY AND MODEL 
A. Preview 


In the following we suggest a way to approximate the T-time in attosecond experiment based on simple mathematical and 
quantum mechanical rules. Our start is a model of Augst et al. [14, 15], where the appearance intensity of a laser pulse for the 
ionization of the noble gases is predicted. The appearance intensity is defined [14] as the intensity at which a small number of 
ions is produced. In this model (in atomic units) the effective potential of the atom-laser system is given by 


V^ffix) = V(x)-xF = -^ -xF, 


( 2 ) 


where F = Fm is the field strength at maximum of the laser pulse (in this work in all our formulas F stands for Fm)-, and 
Zeff is the effective nuclear charge that can be found by treating the (active) electron orbital as hydrogen-like, similarly to 
the well-known single-active-electron (SAE) model [16, 17]. The choice of Zeff is easy recognized for many electron system 
and well-known in atomic, molecular and plasma physics [18-20]. We take one dimensional model along the x-axis as justified 
by Klaiber and Yakaboylu et al [22, 23]. Augst et al [14] calculated the position of the barrier maximum Xm by setting 
dVeff{x)/dx = 0 ^ Xa = Xm{Fa) = {yjZ^ff jF q) and by equating Veff{xm) to the ionization potential Ve//(^a) = —Ip 
(compare fig 1, the lower green curve) they found an expression for the atomic field strength Fa, 




XaF — I'D Fa — 


72 

P 


(3) 


Xa 

and the appearance intensity la = F^. Now we take this idea and relate our argumentation to this model for F < Fa, i.e. for 
the tunnel ionization in the regime of the well-known Keldysh [24] parameter ■yx <1- It is easy to see (see fig 1) that the 
tunnel exit obeys Xe,+ (7^) > XmiF) (note the equality is valid for F = Fa, and for the subscript -h see below) and the energy 
of the tunneling electron is not sufficient to reach the the top of the barrier (as for Fa) suggesting an energy uncertainty, which 
is determined by the energy the electron needs to overcome the barrier and appears in the continuum at the exit point. It 
appears with zero velocity at the exit point Xe,+ according to the strong-held approximation (SEA) due to Keldysh-Faisal-Reiss 
theory [24-26]. Indeed the barrier height at a position x is, compare hg 1: 


hB{x) =1 hB{x) 1 = 1 E - (x) 1 = 1 -Ip + + xF I 


(4) 


that is equal to the difference between the ionization potential and effective potential Veff{x) of the system (atom+laser) at 
the position x, where F = —Ip is the binding energy of the electron before interacting with the laser. Note we can also get Xm 
and the maximum hB{xm) from the derivative of eq (4), dhsix)/dx = 0. An immediately arising question is, what about its 
maximum hB(_Xm) and the energy uncertainty when the electron passes the barrier and is shifted to the continuum or a “quasi“ 
energy level?. First in regard to the derivation of Xa = X'm{Fa) and the eq (3) of Augst et al, it turns out (compare hg 1) that 
the maximum of the barrier height hB{xm) for arbitrary hied strength lies at Xm{F) = yjZ^ff /F. This follows immediately 
from the fact that Xm is determined by the maximum of the effective potential energy for arbitrary hied strength and that is 
the intersection point of the two potentials, V{x) = — and —xF (see hg 1), then 



-xF ^ Xm{F) 



( 5 ) 



4 



Figure 1. 

exit“ point Xe 


Graphic display of the potential curves, the barrier width and the two inner and outer points = (Ip ± 5z)l'2F, the "classical 

c — IplF and the X'm(F) — ^ Z ^ ^ / F position at maximum of the barrier height, {note x^ — ^m(F — F^)) see text. 


Otherwise (to the both sides) one of the two potentials or — xF) slopes down stronger than the other slops up (which 

can be easily gathered from fig 1 and eq (5)), leading to ( 2 ^) < ^e// i^rn) for x 7 ^ Xm- Indeed eq (3) can be generalized as 
the following, for a filed strength F < Fa we get 


F < 


P 

p 


4Z. 


e// 


^ = Ip - 4-^e//^ > 0 


( 6 ) 


The equality ^2 = 0 is valid for F = Fa- We will see that 5^ = ^zi^) = y^p ~ is a key quantity, it controls the 

tunneling process, and determines the time "delay” under the barrier Tj. ^ and the total or the symmetrical T-time , 

subsec. lie. From fig 1 we see that the barrier height at Xm‘ 


hB(Xm) =\-Ip - Vi,ff{x,n)\ = \-Ip + ^j4:Z^ffF\ (7) 

This is the maximum of the barrier height, and by setting hB(xm) = 0 one obtains Fa = 7p/(4Zeyy), which is equivalent to the 
setting Veff{xm) = —Ip as done by Augst et el [14], and can be easily verified from eqs (4), (7). Here we indicate one of the 
failure of the Keldysh time primary resulting form its inadequate definition. If we recall the definition of the Keldysh time, the 
time it takes a classical electron (with an average velocity) to cross a static barrier of a length I [27]. For hs = 0 I = 0 we get 
the T-time = 0 (meaning the ejection of an electron happen instantaneously at F = Fa) because the barrier width vanishes 
for F ^ Fa and (re^_p-£Ce,—1 = 0. But we know, at appearance intensity (7a = T^a) ionization time is equal to ^ (in au) [28] 
(chap. 8) and is not zero. Which seems naturally because the energy gap that has been overcome is 7p, and as we will see later 
this follows immediately from our model, r^^dil^a) = a^nd to a total time T'p^symil^a) = -^1 sae below. As a consequence 
Keldysh time represents a laboratory clock (parametric time), whereas in our following model and time-relation(s) the time 
is dynamically connected to the system (to observe the time form within the system and consider the quantum nature of the 
particle), thus it represents a quantum clock. 


B. Sketch of the model 

The (tunneling-) ionization happens according to SFA with zero kinetic energy at the exit point Xe- Our idea is that the 
uncertainty in the energy can be quantitatively discerned from the atomic potential energy at the exit point AE ~| V{xe.) |=| 

I for arbitrary field strength F < Fa. Then when the electron moves under the barrier in the a^-direction [22, 23], 
its kinetic energy getting smaller, at the same time it moves upwards on the potential energy scale losing potential energy 
(—getting smaller in absolute value). The change happens simultaneously in the potential and the kinetic energy while 
staying at the level —7p-line when tunneling, compare fig 1 , until its kinetic energy becomes zero at the exit point, although 
its (atomic) potential energy not — 7 ^ 0. This can be also gathered from the analysis of the short-range Yukawa and 
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long-range Coulomb potentials given by Torlina et al. [29]. Their conclusion supported with ab initio numerical tests, is that 
for long range potentials, ionization is not yet completed at the ”moment“ the electron exits the tunneling barrier in contrast 
to the usual assumption that ionization is completed once the electron emerges from the barrier. Indeed it is not difficult to 
see that the electric field of the laser pulse shift the electrons along the x-axis direction [22, 23], compare fig 1, reaching the 
exit point with zero velocity, i.e. the electron is forced by the electric field to take and move along a preferred direction, and to 
reduce its kinetic energy to zero at th exit point Xe (the field interacts only kinematic-ally with the electron since no photon 
absorbing), where it is still underling the attraction of the atomic potential V{xe) = — , that defines the uncertainty in 

the energy and acts as a shutter open/closed like in the photon box Gedanken experiment with an uncertainty proportional to 
/\.E '*^1 V{xe) |. As we will see in IV the result is very convincing. 

First it is straightforward to show from eq (4) that for the classical exit point Xe,c^ setting AE ~| V(a^e,c) | leads to incorrect 
T-time (i.e. it fails to predict T-time measured by the experiment). This is because the atomic potential is neglected in 
calculating Xe,c [30]: 

v1 - Vq = 0 - 2Ip = -2F{Xe - xo) Xe,cF = Ip 


where Xe ~ (a^e — a^o), a:o ~ 0 is the initial point of the electron and assuming the electron moves along the x-axis direction 
[22, 23]. It is easy to see, from fig 1, that the ’’classical exit“ point is the intersection of the electric field line —xF with the 
ionization potential of the electron —7p-line, hence Xe^c = which is far from the “correct" exit point. Therefor to use the 
’’classical exit" point, Xe,c-, to determine the T-time will never give a correct answer. It is important to use a correct exit point. 
As seen in fig 1 the relation to other intersection points is simple. The crossing points Xe,± of V^^j^{x) with the —7p-line are 
given by hB{x) = 0, which leads to 


^e,± 


Ip di 5z 
2F 


X 




= Xe,c - Xe,- 


( 8 ) 


where 5z = ^z(E) is given in eq (6), we emphasize the dependence of 5z on Z^ff as done by [20]. Note the origin of the axes is 
at 0. We will see later that using the ’’classical exit” point Xe,c in the uncertainty relation gives the first term in the expansion 
of the T-time obtained using the ’’correct” points Xe,±, eq (8). 

A first arising question in our model is, what happens in the limit of the appearance intensity, i.e. for F ^ Fa = Ip I (AZeff), 
where the electron is shifted from the ground state Eq = —Ip to Ef ^ 0 appearing at Xe,— = Xa = Xe,+ with zero velocity. 
Its energy uncertainty (since no photon absorption) is then (according to our model) A_E'(7^a) ^\ — Zeff/xa\= One sees 
that for atomic field strength (Fa) the electron is heavily disturbed but appearing not far from the nucleus at Xa = XmiEa) = 
with an ionization time that follows immediately from TEUR in eq (1): 



_ _ 1 _ 1 
“ 2 |A£;(Fa)| “ Tp 

as it should be for the ionization process at the Fa [28] (chap. 8). However we will see later that Fa is a special case because 
Xe,— = Xa = Xe,+ (a double solution of h^ix) = 0) and the limit rx^symiF Fa) = ^ is a sum of two equally terms (see 
sec. lie, and III 1). In conclusion, our model is meaningful and the atomic potential energy of the electron at the tunnel exit 
Xe (instead of the gravitational potential in the Einstein-Bohr Gedanken experiment) amounts to calculate the uncertainty of 
the energy in the tunneling process, and hence the T-time by the virtue of eq (1) which leads to an excellent result as we will 
see sec. IV. 


C. Tunneling time 

Our goal now to find an expression to calculate the T-time, and what we need is the correct exit point, where many 
approximations exist. The most used one in the literature can be obtained from eq (4), where the barrier height at the entrance 
and the exit vanishes hB(xe,±) = 0, they are the crossing points of Veff curves with the —7p-line, see fig 1, 

^eff 


hB(x) = 0 


X 


- xF = -Ip 


( 10 ) 
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Solving eq (10) gives immediately eq (8). As seen in fig 1, Xe,— is the inner crossing point, the “entrance” point, and Xe,+ is 
the outer crossing point, the “exit“ point. Physically is argued that the electron escapes the barrier at when it moves in 

the direction a^e,— ^e,+ ^-nd vice versus for the opposite direction, and we will see this presents a useful symmetry property 

of the tunneling process when deriving the T-time. 

Now we can calculate the uncertainty in the energy /S.E[xe,+) by using the exit point Xe,-\- (in the direction Xe,- ^e,+) 
and from this the T-time. From eqs (8), (1) and according to our model we get: 


AE{x,,+) = \^^\= ~ 


Xe,+ (Ip+&z) 

1 1 1 


,unsy — 


(11) 

( 12 ) 


2 AE {Ip — Sz) 

which we call the unsymmetrical T-time rj'^unsy We will see later that a factor (1/2) is missing that can be recovered by 
a symmetry consideration. We first show that the first order of eq (12) is equal to the T-time resulting from using a;e,c, the 
’’classical exit” point, and then we look to the symmetry of the tunneling process. Expanding eq (12) in term of ^ = (^), we 
get in the first order Trp the T-time at the “classical exit” point Xe^c as mentioned above, then 




) = Ie = 

2F 2 AEc 


(13) 


where AEc =l 5 p^l- What about the inner point Xc^ — 1 we could assume, due to the (S^-symmetry between Xc,+ , Xc,-, that 
the electron enter the barrier backwards from Xe,-\- ”entrance“ to Xe,— ’’exit” with an uncertainty (according to our model) 
AE{xe, — ) =\ which leads to (compare eq (11), (12)) 


-1- ^ Cfp—5^ ^ -1- 

2A_E(a;e,-) AZ^ff-F {Ip + Sz) 

This symmetry is deduced in a similar way the Aharanov-Bohm time operator [31] is defined for a free particle T = ^{xp~^ -h 
p~^x) or in more elaborate (the so-called bilinear form) and detailed treatment given by Olkhovsky and Recami [32, 33]. Such 
operators (given by Aharanov-Bohm or Olkhovsky) have the property of maximally symmetric in the case of the continuous 
energy spectra, and the property of quasi-self-adjoint [34] operators in the case of the discrete energy spectra [32] (and the 
references therein), and are the nearest best thing to self-adjoined operators and satisfy the conjugate relation with the 
Hamiltonian and therefore implies an ordinary TEUR [32] and [3] (chap. 1). We use this property, i.e. we assume that the 
maximally symmetric (or almost self-adjoint [33]) property holds for the T-time, for more details see [32], which leads using 
eqs (12), (14) to a simple relation for what we call the symmetrical T-time given by: 


,sym — H” '^p^^ — '^T,i '^T,d 

2 (AB- ^ AE+) 

1 ^ ^ Ip 

{Ip-5c)) ^ZcffF 

where we defined ^ = 1/(2AE'^) = (2(/p di 5z))~^, or (1/2)A_E'='= = AE{xe,±). We call the first term of eq (15) ^ (= ^_) 

and the second term ^p d(~ '^t - ) ^ reason that will be clear in the next section, sec. III. Relation (15) has again (clearly 

because 6z = 0) the correct limit for atomic field strength (compare eq (9) and the discussion after it): 


- ( —^— 

2 \ ilp + 5c) 


'^p ,sys{E —>■ Ea) 


1 1 
Zfp 27/) 


4 4^77 


1 

Tp 


Note that the limit F ^ Fa gives Xe,+ = ^e,— = ^m(^a) = which means that the two points coincide at the top of the 
barrier (a double solution of eq 10). The question is whether this means a symmetry break of the tunneling process? because 
at this limit there is only one direction, namely towards the continuum, so that the “tunneling“ becomes a ’’real” ionization (or 
an ejection) process at the appearance intensity la = F^, and Sz becomes imaginary for super-atomic field strength F > Fa 
(whereas F < Fa is called the subatomic field strength), see below sec IV. A reasonable question in this case, whether the 
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time to reach the “entrance“ of the barrier and under the barrier, intrinsic or dynamical to be measured by a quantum clock, 
becomes after the tunneling a classical, external (or parametric) time due to the break of some symmetry property? so that 
(only) under such a symmetry break a quantum clock (the internal time) coincides with a laboratory clock (the external time). 

In sec. IV we will see that eq (15) (especially gives an excellent agreement with the experimental result of [1]. It 

worthwhile to mention that the T-time measured by the experiment is the time that the electron spends when it moves under 
the barrier from Xe,- to Xe,+, which corresponds to the second term of eq (15) ^ ^ (ip —5 ) ’ further sec. IV result and 

discussion. 


III. PHYSICAL REASONING 

1. Tunneling time and a model of a shutter 

The theoretical mathematical model developed in sec. II can be supported and derived by physical arguments, the only 
difference is that we try to figure a physical insight that helps us to put physics in mathematical relations. From fig 1 and eq 
(8) we get (dg(F) is the barrier width): 


dB{F) = Xe,+ - Xe,- = 


(Ip+ 5,) (Ip-5,) 


(16) 


2F 2F F 

from this it follows 5^ = {xe,+ — a;e, —) F. Because for atomic field strength F = Fa => = 0, dg = 0, we can interpret 5^ 

as the field’s (kinetic) energy exerting in the tunneling process between the “entrance" (or the “inner") point and the ’’outer" 
exit point. The uncertainty in the energy as the electron moves on the —7p-level to a:e,+, i.e. when tunneling the barrier and 
escapes at the exit to the continuum is then Ai?"*" = abs{—Ip + 5^). That means the barrier itself causes a delaying time 
relative to the atomic field strength Fa (5z = 0). The T-time is then obtained from eq (1), 

1 1 


2AE+ 2(/p-<5z) 


(for F < Fa) 


(17) 


which we call ^ (compare eq 15) meaning that is the time duration (time interval) to cross the distance under the barrier 
(in the direction x^^— x^ and escapes at the exit point Xe,+ to the continuum. The term rT,d(F —>■ Fa) = at the 
limit F = Fa account for the turning off the wave packet (or shack-off step, see discussion below) at the ’’entrance-exit” point 
Xa to escape to the continuum, which indicates that the shack-off step or turning off at Xe,+ for F < Fa happens with a time 
delay as given in eq (17). Expanding Sz = /p->/l — {AZ^ffF/Ip‘^) and taking the first we get eq (15), r 4ZeffF ~ sym- 
It means using the symmetrization gives a linearized time duration (linearized T-time). As we already mentioned the limit of 
eq (17), l/(27p), is only one part of the total tunneling time l//p. 

Now we argue that, this picture fits well in the Gedanken experiment of Einstein and the (intrinsic) time r^, ^ (the second 
term in eq (15) or (18)) is rather the time of passage with the shutter open/closed time interval (generated in the internal 
time [5]) is related by the virtue of eq (1) to the uncertainty (1/2)AE'*" = AE = V(xe,+) (note we recovered in eq (15), (17) 
the factor (1/2) missing in eq (12)), and we think that the attosecond experiment with the help of our model represents a 
realization of the photon-box Gedanken experiment (with the electron as a particle instead of the photon) with an uncertainty 
being determined from the (Coulomb) atomic potential due to the disturbing by the field F, instead of (the disturbing by) the 
wighting process and as a result an uncertainty in the gravitational potential [5], as done by the famous prove of Bohr to the 
uncertainty (or indeterminacy) of time in the photon-box Gedanken experiment [5—7]. 


2. Totale time. (F = Fa) 

At the moment one can obtain the total time, i.e. including the time to reach the entrance of the barrier Xe, — , by adding 
the term l/(27p) 

1 / 1 1 \ 

( 18 ) 


2 \lp abs{-Ip + (5z) 

1/1 1 



where the index t only to distinguish between different notations. Nevertheless the term l/(27p) we added is exact only when 
F = Fa or Xe,— {Fa) = Xa = Xe,+ (-Fa), hence we call it r(xi,Xa) 

T{xi,Xa)=^ (19) 

Zip 

where Xi is the initial point. It can be viewed as the response time of the electron to the field, that is, the electron received 
a kick by the field, and is polarized along the field direction, and (jumping up [30]) is moving from xi to the “entrance-exit“ 
point Xa to the continuum, Xi Xg, —(= Xa = Xe^+) this case {h^^F = Fa) = 0, (Sz = 0 ) the most probable 

’’tunneling” mechanism is that the laser field distorts the electron, shakes it up (moving from Xi to Xa), and shakes it off 
(moving to the continuum) at a (total) time given in eqs (15) or (18) ^ = ^. In this model, for F = Fa the 

(illustrative) two steps are not strictly separated, whereas, see further below, for F < Fa they are well separated due to the 
barrier dsiF) > 0 = ^^(-^a)- 


3. Total time for subatomic field F < Fa 

However, Xa = Xe,— = Xe,+ is the maximum “entrance^ point (most right see fig 1), the electron is less disturbed for F < Fa 
and moved to a point Xe,— < Xa that is closer to the initial point Xi [35], this shortens the time to reach the “entrance^ point, 
and we expect that the response of the electron to a small field strength F is weaker than to a stronger field F ^ Fa- The 
time reduction in r^xi, ^e,-) for F comparing to T(xi^Xa) = 1/(2/p) for Fa-, eq (19), is a factor depending on 5z (see discussion 
after eq (16)), because the kinetic energy experiences a change proportional to (x^ _ — xi) F < {xa — xi) Fa- A (weaker) field 
F < Fa not sufficient to compensate for the kinetic energy at the (shake-up) step, instead the electron is at the ’’entrance” 
Xe,— < Xa with a velocity that is sufficient to inter under the barrier and reaches the exit point Xe,+ with zero velocity. Here 
is again an another failure of the Keldysh time (see also discussion after eq (7)) that is the electron did not inter the barrier 
with its initial velocity y^27p suggested by many authors [30]. Now we give the following relation for the time needed to reach 
the ”entrance“ point Xe, — -, and show an explanation further below. 


TT,i = 


1 


1 




2{IpF5z) 2AE- 


( 20 ) 


The factor (comparing to eq (19) for Fa) in the denominator results from two parts. Indeed we follow [5] in that the 
uncertainty is a result of the different reactions or responses of the electron to different filed strengths. The one part comes 
from the difference of moving along the x—axis. i.e. the difference in shifting the electron to Xa with Fa or to Xe,- < Xa 
for F < Fa, which leads to Ai = XaFa — Xe, — F = ^ ~ other part can be deduced from the change 

on the vertical (potential energy) scale. When the electron receives a kick, changing its potential (on a vertical scale), its 
atomic potential has a different changes between V(xi) and V(x) for x = Xa or Xe, — - But this is a result of different responses 
(on energy scale), while the electron is forced to follow an orientation along the (opposite) field direction at Xa, or at Xe, — . 
Therefor this part can be approximated from the difference AK(a;) in the atomic potential at the different “entrance‘ points, 
which gives A 2 = V(xa) — V{xe, — ) = ^ (Jp+'^z) _ ^ ^ difference equal to 

Ai -h A 2 = Sz between reaching the ’’entrance” point Xe,— (for F) relative to the “entrance-exit“ point Xa (for Fa), leading to 
an energy uncertainty AE = abs{—Ip — 6z) = (/& + 5), for F < Fa, and with a time rj'^i{F) = ^ for arbitrary subatomic 

field strength F, hence one obtains eq (20) instead of eq. (19). We have explained so far eq (17) and eq (20), in doing so 
we explained the physical meaning of the symmetry consideration done above (see eq (15) and the discussion before) and we 
obtain from eqs (17), (20) the result obtained in eq (15): 


= '^T,i + '^T,d = '^T,+ + '^T,- (21) 

1 ^ _ 4 

{Ip-Sz)J AZeffF 

the first term = 1/(2E~) corresponds to the first step, where the electron is shacked up and moved to the ”entrance“ Xe,- 
(or Xa for Fa) that takes the times Tj.. = Trp The second term _ = 1/(2F'*") corresponds to the actual T-time or moving 


2 \dp 
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under the barrier and shaken off to the continuum, with a time delaying ^ _ due to the barrier relative to the atomic 

field strength Fa , where 5z = ^ and rj. = . 

For Fa, as already mentioned, the second or shake-off step follows immediately the first or shack-up step and the two steps 
are not strictly separated. For F < Fa the two steps model of the tunneling process are well separated. They happen with 
opposite directions at the time scale, the first step is less time consuming since F causes a smaller disturbance relative to Fa, 
and the the electron moves not far from its initial position for small F, Xe,— < Xa, whereas the second step happens with a 
time delaying, Xe,+ > Xa, relative to the ionization at atomic field strength. So far our theoretical model is assisted with an 
explanation through a physical reasoning. In the following we show and discuss our result for He-atom with a comparison to 
the experiment [1, 2]. 


IV. RESULT AND DISCUSSION 

In fig 2 we show the results of eq (12) the unsymmetrical ansy ’ symmetrical T-time Tj. . The results 

for Tj, ^ eq (17), and again the the symmetrical (or total) T-time of eq (21) are shown in fig 3. Note eq (21) is identical 

with eq (15), whereas eq (17), is the second term of eq (21) or (15), which is the actual T-time, i.e. the time needed to pass the 
under barrier region (a^e,- ^e,+) Q-nd escape at the exit point to the continuum, and that is usually the T-time measured in 

the experiment. The results are for the He-atom in a comparison with the experimental result of [1]. The experimental data 
and the error bars in the figure were kindly sent by A. Landsman [36]. We plotted the relations (12), (15), (17) and (21) at 
the values of the field strength at maximum of the elliptically polarized laser pulse (A = 735, elliptical parameter e = 0.87, 
F = Fo/'s/l + £^) used by the experiment exactly as given in [36]. In fig 2, the upper two curves are the unsymmetrical T-time 
Tt unsy'> (1^)’ different models of the effective nuclear charge Z^ff, that the tunneling electron experiences during 

the tunneling process. Accordingly the lower two curves are the symmetrical T-time, Xj, sym^ different 

models of Z^ff- We mention that eq (18) (not plotted) gives a closer result to eqs (15), (21). The two different effective charge 
models are from Kullie [37], with Z^^j^ ^ = 1.375 and dementi [38] with Z^^j^ ^ = 1.6875 . In fig 2 we see that our unsy 
not close to the experimental data of [1]. Whereas is close for both models of the Zeff- But we notice, see discussion 

further below, that the first term in eqs (15), (21) is much smaller than the second . < Tj. ^ for small F (relative to Fa). 

Concerning Z^ff we see for small field strength F < 0.05 that with ^ closer to the experimental data 

(and especially for see discussion below fig 3). The reason is that the ^ r^odel is a H-atom-like model, based on 

the assumption that the first electron of the He-atom occupies the Is-orbital (with probability density |^(r)p, which screens 
the the nuclear charge and the second electron is treated as an active electron or a “valence“ electron [37] similar to the SAE 
approximation. This is a good approximation when the tunneling electron moves far from the left atomic core (He"^), or the 
barrier width is large a:e,+ > 15 au, hence the better agreement, and possibly this is important for smaller field strength in the 
region where 7 ^ ~ !• In the range of larger field strength multielectron effects are expected and the model -^^yy c based on 
Hartree-Fock calculation is more reliable, where the electron moves not far from the left atomic core (small barrier width) and 
hence the better agreement in this region. It is likely that a model depending on the x-coordinate Zeff{xe,+) will achieve a 
better agreement that smoothly fits the two regions. 

Now we look to fig 3, where eq (21) and ^ eq (17) are shown. Eq (21) is the same as eq (15) (shown in fig 2). For 

the Tj, ^ we see an excellent agreement with the experiment. As already mentioned Trj, ^ corresponds to the T-time measured in 
the experiment, that is the time (interval) needed to move under the barrier from the entrance to the exit point and escape to 
the continuum with a shack-off, or between the instant of orientation at Xe,- an the instant of ionization at Xe,+, which is the 
time spent in the classically forbidden region [1]. For small F < 0.055 au, Z^ff^K gives better agreement with the experiment, 
whereas for larger field strength Z^ff^c is more reliable, where multielectron effects are expected due to the decreasing width 
of the barrier and the tunneling electron is closer to the first one when it traverses the barrier. 

In fig 3 we see that the difference between the total or symmetrical T-time Tj. and the (actual) T-time Tj. ^ is small, 
because the second term ^ in eq (21), incorporates the delaying time caused by the barrier and is the main time of the 
tunneling process for large barrier. Whereas the first part term ^, is due to the shake-up of the electron by the field moving 
it to the ’’entrance” Xg, —, which is small for small F. For large field strength the two parts become closer because the barrier 
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Field strength (au) 

Figure 2. T-time unsy (1^): 3'rid gym (1^)’ two models [37] and [38]. Time is in attosecond units vs laser field strength 

in atomic units, corresponds to the tunneling ionization of the He-Atom in strong field, in good agreement with experimental result [1, 2, 36]. 
Experimental values are kindly sent by A. Landsman [36]. 



0,04 0,05 0,06 0,07 0,0S 0,09 0,1 0,11 

Field strength (au) 


Figure 3. T-time ^ eq (17) and sym (^t)> (21) is identical to eq (15), see fig 2), for two different Z^ff models as in fig 

2. Time vs laser field strength as in fig 2, corresponds to the tunneling ionization of the He-Atom in strong field, in excellent agreement with 
experimental result [1, 2, 36]. Experimental values as in fig 2. 


width getting smaller SzjF = (xe,+ — Xe, — ) 0 and for appearance intensity (Sz = 0) they become equal. Concerning Zeff 

in fig 3, we readily see for Tj, ^ the same behavior as in fig 2. 

In fig 4 we plotted our result ^ eq(17) together with the FPI result of [36] (data were kindly sent by A. Landsman and 
C. Hofmann [36, 39]). From the figure we see a good agreement between the two results, and the difference is smaller than the 
experimental error bars. Indeed we expect that the FPI would agree better for large field strength F > 0.055 with the lower 
curve {Ze.ff,Ci green). For small field strength, FPI is more or less close to both curves (green, blue), but our upper curve 
{Zeff^ki blue) tends to be in a better agreement with the experimental values for smaller field strength. An important point is 
that our model and result(s) predict a (real) T-time (time of passage) of a single particle, it is not distributive of an ensemble 
(although indeterminately in regard to the uncertainty relation) and we make no assumption about the path of the particle 
inside the barrier, whereas the FPI treatment is probabilistic/distributive that makes a use of all possible (classical) paths inside 
the barrier that have a traversal tunneling time r. Furthermore, Landsman et al [40] uses the time tq, which is determined by 
the measurement to coarse-grained the probability distribution of the T-times to achieve the aimed results, although Sokolovski 
[13], [3] (chap. 7) claims (in regard to his FPI description) that no real time is associated with the tunneling. We think that the 
two views are rather complementary as it is usual in quantum mechanics: wave/particle or individual (single particle)/statistical 
(distributive) etc. The Larmer Clock should in principle agrees with our result, the result of Landsman [36] shows that the 
agreement is good only for F ^ 0.06 — 0.1, hence Larmer Clock values are inferior for F < 0.06. The same holds for F > 0.1, 
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Field strength (au) 


Figure 4. T-tinie ^ eq (17), for two different models as in fig 2 together with FPI result [36] and the experimental result [1, 2], 

Experimental values and FPI result are kindly sent by A. Landsman and C. Hoffman [36, 39]. 


although the difference is smaller than the error bars, but the trends of of Larmer Clock curve for large field strength seems 
not correct. In general the Larmer Clock curve is fiat comparing with the other curves and the experimental data, which is 
difficult to understand. 

In fig 5 we show the tunneling time versus the barrier width dsiF). The T-time shows a linear dependence with the 
barrier width dsiF) in the region F = 0.04 — O.llau and a limit l/(2/p) at ds^Fa) = 0. The other limit for very large 
barrier width {F ^ 0, (53 ^ Ip) is 4 z^^fF ~ '^T.sym ~ which is straightforward because for very large barrier 

7'T,d >> '^T,d ~ '^T,sym- We note as Seen in fig 5 that the time spent by a particle (a photon) to traverse the same 

barrier width with the speed of light is much smaller than the T-time of the electron in the He-atom. 

At the limit F = Fa of the sub-atomic field strength the tunneling process is out and an ionization process called the 
’’above the barrier decay“ is beginning. For supper-atomic field strength F > Fa-, Sz becomes imaginary (and so the crossing 
points, compare eq (8), but still a real Xm = -s/Zeff/F), which indicates that the real part of ^ or is the limit 

for ’’real” time tunnel-ionization process. Indeed in this case the atomic potential is heavily disturbed and the imaginary part 
of the time ^ is then due to the release of or escaping the electron (at Xm{F)) from a lower energy level than —Ip (and 
possibly escaping with a high velocity), where the ionization happens mainly by a shack-off step [28] (chap. 9). Here we see 
the clear difference between the quantum mechanical and the classical clocks [5, 31]. Classically we can make the interaction 
time with the system arbitrarily small, the real part of the time can be made arbitrary small, and an imaginary part is absent. 
In quantum mechanical the tunneling-ionization time has a real part limit ^ = l/(27p), an imaginary part arise when the 
field strength is larger than the atomic field strength Fa , in both terms . and ^. 

However in our treatment, although, ^ have an imaginary part, we get a real total or symmetrical T-time Trp^sym = 

ionization processes with an arbitrary field strength. It becomes very small for very large field strength, and probably 
it loses its validity in this regime suspecting a break of some symmetry and non-linear effects arise, and the interaction becomes 
physically a different character. It certainly also loses its validity in the multiphoton regime, i.e. for large Keldysh parameter 
7 >> 1, where F << Fa- It is apparent from Tp^sym, eqs (15), (21) that the T-time has no imaginary part when the symmetry 
of the time is considered, i.e. when assuming the maximally symmetrical (or quasi-self-adjoint) property discussed in details by 
Olkhovsky [32]. It is now the question to which extend the above relation ^z~^jF‘’ Preserves its validity for F > Fa (or for small 
F « Fa, where ‘Jk « 1), where or what is/are the limit(s) of its validity? Probably a break of some symmetry for F > Fa 
(or F << Fa) can give a hint to answer this question. Finally we mention that for F > Fa or intensities I > la Stark-shift, 
relativistic and non-linear effects become large, the perturbation theory breaks down (which is valid for small parameter FjFa 
[28]) and several regions appear at intensities larger than the appearance intensity la, such as the critical Ic and the saturation 
Is intensities. Is > Ic > la, where multiple ionization occurs [28] (chap. 7, 9). 

Conclusion We presented in this work an analysis for the tunneling process and the tunneling time in attosecond 
experiment and found an accurate and simple relation to calculate the tunneling time for the important case of He-atom, where 
a reliable experimental data are available. Our result (especially the T-time r^^) was shown to be in excellent agreement 
with the experiment [1, 2] and with the Feynman path integral treatment of [36] although for small field strength our result of 
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Barrier width (au) 

Figure 5. T-time Tj, ^ eq (17) in as vs barrier width dQ{F) in au. for two different models as in fig 2. The lines at the bottom of the 

figure show the time spent by a particle (a photon) traversing the same barrier with the speed of light. 


Zeff.K tends to agree better with experimental values. Note in figs 2-5 we use for the evaluation of our result the same values 
of the field strengths used by the experiment, i.e the filed strength at maximum, see [36]. The T-time in our treatment is 
dynamical or intrinsic-time type and represents a quantum clock, i.e. to observe the time form within the system and consider 
the quantum nature of the (bound) particle, in contrast to the classical Keldysh time which is external or parametric, where 
we indicated (one or two of) its failures to treat the T-time in our (study) case. 

Further we suggest a model of a shutter to the tunneling process in attosecond experiment, and we think the experiment 
together with our tunneling model (subsec. IIB, IIC) offer a realization of the Bohr-Einstein’s photon box Gedanken experiment, 
with the electron as a particle instead of the photon and with the uncertainty being determined from the (Coulomb) atomic 
potential instead of the gravitational potential. Our treatment suggests that a symmetry (maximally symmetrical or quasi- 
self-adjoint) [32] assumption to calculate the T-time is important and gives a hint to the search for a time operator in the 
tunneling process and maybe for a general time operator in quantum mechanics. Our result uses two models of the effective 
charge Zeff of the left core He'^ that the tunneling electron experiences. The Z^^^ ^ = 1.375 of Kullie [37], that based on a 
similar model to the single-active-electron, is better for small field strength F < 0.055, whereas Z^^^ ^ = 1.6875 of dementi et 
al [38] is more reliable for larger field strength, because it is based on the Hartree-Fock calculation, and that is justified, when 
the multielectron effects are not negligible. 
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